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An iterative solution best suited for a Monte Carlo implementation is presented for the non-forward 
BFKL equation in a generic color representation. We introduce running coupling effects compatible 
with bootstrap to all orders in perturbation theory. A numerical analysis is given showing a smooth 
transition from a hard to a soft pomeron when accounting for running effects. 



1 Introduction 

The solution to the LL BFKL equation projected in the color singlet in the non-forward case was first 
calculated in the seventies by Balitski, Fadin, Kuraev and Lipatov |lj|5}. Three decades later a solution 
based on a sum of iterations of the kernel in transverse momentum space was proposed at LL (6| and 
Next to Leading Logarithmic(NLL) ||7) accuracy and used to build up a Monte Carlo code to provide 
numerical studies of it. Similar studies for the color octet representation were given in J8- 10 1. 



The present work is based on the formalism given in [6|. The study is extended to give a solution 
in a general color group representation 1Z and presented in a way such that all the infrared divergences 
appear as an overall factor in the gluon Green function so that a finite piece can be identified and treated 
numerically. In order to account for higher order corrections and study the properties of the infrared the 
running of the coupling is introduced. How to account for it is a non well-defined problem. There is no 
theoretical strong restriction on it and different possibilities have been suggested in the literature fTT| - 
[16). We insert it in a way consistent with gluon reggeization, as proposed in 1 14 -T6J, which naturally 



leads to the appearance of renormalon singularities in the infrared. The fact that the solution is given 
in transverse momentum and rapidity space makes possible to study diffusion properties, analyzed in 
Sec. 1 3 - 1 1 The aim of this review is to give a very short qualitative explanation of the main points of the 
work presented. We refer the reader to the article in preparation [ 17 1 and the references given here for 
the calculation. 



2 Non-forward BFKL equation in a generic color representation 

The infrared divergences that appear in a general color representation can be written as an overall factor 
in the gluon Green function. To show this one has to regularize half of the divergences in the gluon 
Regge trajectory using dimensional regularization (D = 4 — 2 e) and introduce a mass parameter A 
for the remaining ones. In doing it, the dependence on A cancels out with the contribution of the real 
emissions while the one on e remains in the factorized term, leading to a solution to the non-forward 
BFKL equation independent of A for A — » 0. 
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The divergent term depends on the choice of the color projection. In the singlet representation, for 
example, Cx = 1 and all infrared divergences cancel out. The non-forward equation for the color singlet 
can be solved in an iterative way by going back to momentum space using a Mellin transform. Details on 
how this is down are given in [7]. The result obtained is the product of an exponential term depending on 
the A and 1 /e parameters and a finite part that we denote by H (q x , q2 ; q; Y) . We analyze this function 
in sec. [3] It is convenient for the numerical study to express the gluon Green function .F(qi, q2; Y) as 
a function of the azimuthal angle between the two-dimensional vectors qi and q2, its Fourier conjugate 
variable or conformal spin n and the anomalous dimension 7. 



2.1 Infrared effects 

A way to introduce running coupling effects at in the analysis of a 2 — > 2 partonic process at LL accuracy 
is to replace the reggeized gluon propagating in the t-channel by a gluon (or renormalon) chain fl8) . 
Our choice of accounting for the running is based on this approach. A comparison of the LL gluon 
trajectory with the new one tells us that the only needed change to be done in the analytic expressions 
is the replacement k 2 — > T](k), with r/(k) = k 2 /a s (k 2 ). In order to define the new BFKL kernel the 
bootstrap condition is imposed so that gluon reggeization is still justified. This procedure naturally leads 
to the appearance of renormalons as power corrections that could let us learn about the properties of the 
infrared. Details on this setup can be found in |TT4]-fTSI . 

Concerning the choice for the running we use a parametrization which freezes in the infrared and it 
is consistent with global data of infrared power corrections to perturbative observables [ 19 1. When the 
external transverse momentum scales in the gluon Green function are perturbative enough, this model 
for the running cannot be distinguished with a perturbative one with a Landau pole. Nonetheless, we do 
find sensitivity to the IR finite model for sufficient small values of these scales. 



3 Numerical analysis for the color singlet 



Figures 2a show the convergence of the sum defining 



the function % (qi, q2; q; Y). In can be seen how for a 
fixed value of Y and the coupling a s a finite number of 
terms in the sum is enough to have a good accuracy for 
the gluon Green function. As the value of the effective 
parameter a s Y gets larger the Green function is more 
sensitive to high multiplicity terms, following a Pois- 
sonian distribution. It can also be noticed how the dis- 
tribution in the number of iterations of the kernel gets 
broader for larger center of mass energies although the 
convergence is always good. 

It is also instructive to study the solution in terms 
of the different Fourier components in the azimuthal 
angle between the two momenta q! and q 2 . A com- 
plete analysis is shown in Figs. 2b for both forward 
and non-. It can be seen how the only rising component is the 




Figure 1 : Diffusion pattern for fixed coupling 
(blue, dotted line) and running coupling (red, 
solid line). 



2c 



n = one. For completeness, fig 
shows the dependence of the solution on the azimuthal angle for the sum of all Fourier components. The 
collinear limit is investigated in Fig. [2d] 
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3.1 Diffusion 



The diffusion [20. 21 1 of the transverse scales in the BFKL ladder has been studied in terms of the 
average value (r) of r = log ^(qi + X) ki) 2 ^ as a function of the rapidity Y' along the gluon ladder. 
For each set [qi , q2 , Y] (where qi and q2 are the transverse momenta of the edges of the ladder) the 
non-forward BFKL equation is solved numerically allowing the study of the evolution of (r) along the 
ladder as well as the weight of each configuration point in n-momenta phase space to the total solution. 
Fig. [T] compares the diffusion pattern found for a fixed choice of the strong coupling and a version 
with running coupling inserted as explained in the previous section. The straight line in the middle 
corresponds to (r) while the upper/lower curves are the mean plus/minus the standard deviation. The 
set of values taken for the plot are k a — 5 GeV, k\, = 4 GeV and Y = 1. The figure shows how the 
version with running is shifted to the infrared providing a smooth transition from the hard to the soft 
pomeron. 
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Forward LO BFKL Singlet Green function (q = GeV, c R = 1 ) -a Non-forward LO BFKL Singlet Green function (q = 5 GeV, c R = 1 ) 




Number of emissions n -o Number of emissions n 



(a) Distribution in the contributions to the BFKL gluon Green function with a fixed number of iterations of the kernel, 
plotted for different values of the center-of-mass energy, and a fixed a s = 0.2. Left: forward case; right: non-forward 
case with q = 5 GeV. 



Forward LO BFKL Singlet Green function (q = GeV, c R - 1 ) 



Non-forward LO BFKL Singlet Green function (q = 5 GeV, c R = 1) 




n 1 
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(b) Projection of the gluon Green function on different Fourier components in the azimuthal angle between the qi 
and q2 transverse momenta. Left: forward case; right: non-forward case with q = 5 GeV. 




(c) Gluon Green function dependence for the full range in (d) Collinear behavior of the gluon Green function, 

azimuthal angles. 



Figure 2: Numerical analysis 
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